(N 

o 

(N 



in 

(N 



(N 

> 



X 

5-j 



A simple method for the Kramers-Kronig analysis of 
reflectance spectra measured with diamond anvil cell 



H. Okamura 



Department of Physics, Graduate School of Science, Kobe University, Kobe 657-8501, Japan 



E-mail: okamura@kobe-u.ac.jp 



Abstract. When the optical reflectance spectrum of a sample under high pressure is studied 
with a diamond anvil cell, it is measured at a sample/diamond interface. Due to the large 
refractive index of diamond, the resulting reflectance Rd(oj) may substantially differ from that 
C ' measured in vacuum. To obtain optical constants from Rd(uj), therefore, the usual Kramers- 

\ Kronig (KK) analysis cannot be straightforwardly applied, and either a spectral fitting or a 

modified KK transform has been used. Here we describe an alternative method to perform KK 
' analysis on Rd(uj). This method relies on the usual KK transform with an appropriate cutoff 

and extrapolation to Rd(uj), and may offer a simpler approach to obtain infrared conductivity 
from measured Rd(uj). 
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1. Introduction 

Infrared (IR) spectroscopy has been a powerful tool to study the microscopic carrier dynamics 
and electronic structures in strongly correlated electron materials, such as rare earth (f electron), 
transition metal (d electron), and organic (p electron) compounds [I]. The IR spectroscopy 
technique has been also performed under high pressure using a diamond anvil cell (DAC) [2-13] 
since the strongly correlated materials show many interesting physical properties under high 
pressure. In a DAC, a pair of diamond anvils and a thin metal gasket are used to seal a sample 
| and a pressure transmitting medium [6]. A typical diameter of the diamond surface is 0.8 mm to 

reach a pressure of 10 GPa, and 0.6 mm to reach 20 GPa. Therefore the sample in this experiment 
should have dimensions of the order of 100 /um. To perform an infrared (IR) reflectance study on 
such a small sample under the restricted sample space in a DAC, synchrotron radiation (SR) has 
been used as a bright source of both far and mid-infrared. In fact, high pressure IR spectroscopy 
with DAC is currently one of the major applications of IR-SR [4-7,9-13]. 

With a DAC, the reflectance is measured between the sample/diamond interface, in contrast 
to the usual case of sample/vacuum or sample/air interface. The normal- incidence reflectance of 
a sample relative to a transparent medium of (real) refractive index no is given by the Fresnel's 
formula as [TH [15] : 

y J (n + n ) 2 + k 2 w 

Here, h = n + ik is the complex refractive index of the sample, and no=2.4 for diamond and 
1.0 for vacuum. Hereafter, we denote R(uj) at sample/diamond interface as Rd(co), and that 
at sample/vacuum interface as Ro(td). From Eq. (1), it is easily seen that i?d(w) of a sample 
measured in DAC may be substantially different from Rq(uj). 



The purpose of this study is to consider the Kramers-Kronig (KK) analysis of Rd(u) data 
measured in DAC. KK analysis has been widely used to derive optical constants such as the 
refractive index, dielectric function and optical conductivity from a measured Ro(uj) spectrum 
|14U15j . However, due to the difference between Rq(uj) and Rd(w) discussed above, the usual KK 
analysis method cannot be straightforwardly applied to i?d( w ) [IS]- To derive optical constants 
from Rd(u), therefore, previous high pressure IR studies used either a Drude-Lorentz spectral 
fitting [SlElElElElinilElora modified KK transform [3 EJ QUI CE3] • In this work, we propose 
a different method, which relies on the usual KK transform with an appropriate cutoff to the 
Rd(co), as an alternative approach to obtain the infrared <t{lo) from R^{oj). The validity of the 
proposed method is demonstrated with actually measured reflectance data of PrRu4Pi2- 

2. Kramers-Kronig analysis of reflectance spectra 

The complex reflectivity of the electric field, f, is expressed as [14^ [T5] 

fH = ^#4 = ^K' eM - (2) 

v ; n +n(uj) w w 

Here r(oj) is the square root of the reflectance R(co), which is actually measured in experiments. 
Then the real and imaginary parts of n can be expressed in terms of r and 9 as 

1 - r 2 

n = n z • no (3) 

l + r 2 + 2rcos0 u V ' 

and 

— 2rsin# 

fc= l + r 2 + 2rcos^ n °' (4) 

respectively. Therefore, if 0(uj) can be derived from measured r(u>) with KK analysis even for the 
sample/diamond reflection case, n(u) and k{uj) can also be derived simply by setting n$=2A 
in Eqs. (2) and (3). Then, the imaginary part of the complex dielectric function is given as 
62 = 2n(uj)k(uj), and the optical conductivity is given as cr(u) = ^£2(0;) [T4"l PTd]. 
In performing KK analysis on reflectance data, usually the logarithm of f, namely 

lnf(o;) = lnr(w) + i9(cj) (5) 

is regarded as a complex response function. In the case of sample/vacuum reflection, the KK 
relation between lnr(u) and 9(uj) is expressed as |141 115] 

2u [°° InrV) , 
9{u) = P / — zdw . (6) 

7T Jo W — UJ 

Here, P denotes the principal value. In deriving this relation, it is required that lnr(w) has no 
poles in the upper complex Cj plane when \u\ is finite. This is correct since r(cu) — > only when 
oj — >• 00 in the case of n,Q=l. However, when no > 1 as in the case of sample/diamond interface, 
h = no may be satisfied at some point on the upper imaginary axis [TB]. This point is denoted 
asii = i/3, where (3 is a real, positive and finite number. When n = no, f = from Eq. (2) and 
lnr therefore has a pole at Cj = i/3. Accordingly, the KK relation in this case must be modified 

to [E] 

e(oj) = -—P f°° \ dJ + [tt - 2 arctan (B/u)]. (7) 

7T Jo 00 z — 0J A 

Namely, the presence of a medium with no > 1 brings an extra phase shift, indicated by the 
square bracket in Eq. (7), into the KK relation. Note that the extra phase shift is a decreasing 



function of /3/lj, and that the original KK relation of Eq. (6) is recovered when /3/u) —¥ oo 
|16| . Detailed theoretical considerations on the extra phase shift in various situations have been 
reported [T71 [18]. In the case of actual experimental studies, however, the precise value of (3 
may not be known. Accordingly, the value of f3 has been estimated from experimental Rd(w) 
data by use of a combination of DL fitting and the modified KK transform [3j [HJ [TOl [T3] . In this 
method, one uses Eq. (7) with r{oo) = y/R^ito) and looks for a value of (3 that well reproduces 
the ct(uj) given by a DL fitting of R<±(uj)- 

Note that, on the other hand, if the frequency range of interest is lower than the value of f3, 
effects of the extra phase may be only minor, and the usual KK transform of Eq. (6), combined 
with the use of uq=2A in Eqs. (3) and (4), might give sufficiently accurate values of optical 
constants. We will examine the validity of such a procedure in the next section. 

3. Simulation with measured reflectance spectra 

Here, we use Rq{<jo) spectra actually measured on PrRu4Pi2 |19| . This compound is well known 
for showing a metal-to- insulator transition at about 60 K [20], and a clear energy gap in cr(u;) 
was observed in our previous work [19]. Here we use Ro(uj) at 60 K (metal) and 9 K (insulator) 
measured over a wide photon energy range of 0.008-30 eV [19j . which are shown by the blue 
curve in Figs. 1(a) and 2(a). The procedure is the following. 

(i) The full Ro(lu) spectrum is KK analyzed with Eq. (6) to obtain n(u), k{oj) and a(u). 

(ii) The above n(oj) and k(u) are substituted into Eq. (1) with no=2.4 to derive i?d( w ) that is 
expected in a DAC. 

(iii) The Rd(u) obtained above is used with the usual KK transform of Eq. (6) and uq=2A in 
Eqs. (3) and (4), to obtain n(uj), k(u) and a(u). Before this is done, an appropriate cutoff 
and extrapolation are made to the Rd(ui), as described in detail below. 

If the KK analysis on Rd(uj) works properly, the resulting ct(lo) from (iii) should well agree with 
that given by Ro(uj) and the usual KK analysis. 

We first examine the 60 K data. The expected Rd(oj) at 60 K obtained by (i) and (ii) is 
indicated by the black curve in Fig. 1(a). In carrying out the integration in Eq. (6), the Ro(oj) 
spectrum were extrapolated below 0.008 eV and above 30 eV with the Hagen-Rubens and u;~ 4 
functions, respectively. It is seen in Fig. 1(a) that Rd(u) shows very high values above about 
4 eV. This physically unrealistic feature resulted from the unphysical assumption of constant 
no=2.4 in the entire spectral range. In reality, of course, the refractive index of diamond cannot 
be constant and real near and above the band gap, where it shows strong light absorption. In 
addition, when u — > oo, n(u>) — > 1 and k(u) — > 0, and therefore Rd(^) — > must hold just like 
any other material. Accordingly, before performing KK transform in (iii), Rd(u) in Fig. 1(a) 
was cut off at some energy w c , and then it was extrapolated with cj -4 function. Several different 
values of uj c were tried, as shown in Fig. 1(a). For each value of oj c , the usual KK transform 
of Eq. (6) was made to get 0(u), which was then used to derive n(w) and k(u) with no=2.4 in 
Eqs. (3) and (4), and to finally obtain a{ui). 

Figure 1(b) shows the a{ui) spectra obtained as described above, with different values of ui c 
which are also indicated in Fig. 1(a). It is seen that the obtained a{uj) spectra strongly depend 
on uj c . With u> c =2.0 eV, the resulting ct(uj) below 1.5 eV agrees very well with that derived 
from the original Ro(uj). (Actually, w c =2.2 eV gives the best agreement, but w c =2.0 eV data 
is shown instead. This is because the 2.2 eV data almost completely overlaps with that from 
Ro(lo), making it difficult to distinguish them in the figure.) The result for the 9 K data, where 
the sample is an insulator (semiconductor), is also shown in Fig. 2. A good agreement is again 
observed between the <j{uj) derived from the full Ro(lo) and that from Rd(u) with uj c =2 eV. 
The spectral range of their good agreement is below 1.5 eV, which is similar to the case of 60 K 
data discussed above. Note also that both 9 K and 60 K data show good agreement with the 
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Figure 1. (a) Ro(oj) is the reflectance spectrum of PrRu4Pi2 measured at 60 K in vacuum |19| . 
and -Rd( w ) is that expected in a DAC calculated from Ro(uj) as described in the text. The green, 
red, and light blue curves are uj~ a extrapolations with cutoff energies of w c =1.5, 2, and 5 eV, 
respectively, (b) The optical conductivity (a) obtained with KK analysis of Ro(w) is compared 
with those obtained with KK analysis of Rd(^>) with a cutoff at cl> c =1.5, 2, and 5 eV, and u^ 4 
extrapolations above them. Below 1.5 eV, ct{oj) obtained from R^{uj) with u c =2 eV agrees very 
well with that obtained from Rq(uj). 

common value of oj c = 2 eV. These results show that, for any spectral change in Rd(uj) (either 
temperature- or pressure-induced) below 1.5 eV, the corresponding (j{uj) can be obtained by 
the present method. In actual high pressure studies of strongly correlated materials with DAC 
[2-13], i?d(w) is usually measured below 1-2 eV. Hence, above the high energy limit of Rd(u) 
measurement, the Rd(u) expected from Rq(uj) can be connected to the measured Rd(co), with the 
cutoff and extrapolation discussed above. Then the connected -Rd(w) may be KK transformed to 
obtain cr(u), as discussed above. An obvious condition required for this method to work properly 
is that the pressure- and temperature- induced changes of Rd(u) should be limited below certain 
energy, which is 1.5 eV for PrRu4Pi2 as seen in Figs. 1(b) and 2(b). This condition is actually 
met in the high pressure data of PrRu4Pi2, which has enabled us to derive its cr(u) under 
pressure up to 14 GPa using the present method [2Tj . 

We have also done similar simulations for other compounds, both metals and insulators, 
using actually measured data, and have obtained similar results. Namely, when an appropriate 
cutoff and extrapolation are applied to Rd(uj), the usual KK transform of Eq. (6) gave 



1.0 



DC 0.5 



PrRu 4 P 12 



R (co) 
Measured 
in vacuum. 




| 2000 

c 



1 1 1 1 1 1 

. (b) 




i i 

co c =5 eV 






J / Fr ° m R ° (C0) / 








. co c =2 eV K 






co c =1 .5 eV^^/^^-A 




i . . . 






0.01 0.1 

Photon Energy (eV) 



Figure 2. (a) Ro(co) is the reflectance spectrum of PrRii4Pi2 measured at 9 K in vacuum |19j . 
and Rd(w) is that expected in a DAC calculated from Ro(u) as described in the text. The green, 
red, and light blue curves are uj~ a extrapolations with cutoff energies of u> c =1.5, 2, and 5 eV, 
respectively, (b) The optical conductivity (a) obtained with KK analysis of Ro(w) is compared 
with those obtained with KK analysis of Rd(^>) with a cutoff at cl> c =1.5, 2, and 5 eV, and 
extrapolations above them. 



spectra which agreed very well with those directly obtained from the wide range Ro(uj). A 
limitation of the present method is, as already mentioned above, it can give correct <r(u;) only 
below certain photon energy (1.5 eV in the case of PrRu4Pi2). Hence this method is useful 
when the temperature and pressure dependences of Rd(oj) is limited to below certain energy. 
In addition, to use the present method, it is required that i?o( w ) is known over a wide enough 
photon energy range, since n(w) and k(u) must be obtained from Ro(uj) with the usual KK 
analysis. While a mathematically rigorous justification of the proposed method is beyond the 
scope of this work, this method may be very useful as a simple analysis technique of reflectance 
spectra measured under high pressure with DAC. 
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